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ABSTRACT 

RM Synthesis was recently developed as a new tool for the interpretation of polarized 
emission data in order to separate the contributions of different sources lying on the same 
line of sight. Until now the method was mainly applied to discrete sources in Faraday space 
(Faraday screens). Here we consider how to apply RM Synthesis to reconstruct the Faraday 
dispersion function, aiming at the further extraction of information concerning the magnetic 
fields of extended sources, e.g. galaxies. The main attention is given to two related novelties 
in the method, i.e. the symmetry argument in Faraday space and the wavelet technique. 

We give a relation between our method and the previous applications of RM Synthesis to 
point-like sources. We demonstrate that the traditional RM Synthesis for a point-like source 
indirectly implies a symmetry argument and, in this sense, can be considered as a particular 
case of the method presented here. Investigating the applications of RM Synthesis to polar- 
ization details associated with small-scale magnetic fields, we isolate an option which was not 
covered by the ideas of the Burn theory, i.e. using quantities averaged over small-scale fluc- 
tuations of magnetic field and electron density. We describe the contribution of small-scale 
fields in terms of Faraday dispersion and beam depolarization. We consider the complex po- 
larization for RM Synthesis without any averaging over small-scale fluctuations of magnetic 
field and electron density and demonstrate that it allows us to isolate the contribution from 
small-scale field. 

A general conclusion concerning the applicability of RM Synthesis to the interpretation of 
the radio polarization data for extended sources, like galaxies, is that quite severe requirements 
(in particular to the wavelength range covered by observations) are needed to recognize at 
least the principal structure of the Faraday dispersion function. If the wavelength range of 
observations is not adequate we describe which features of this function can be reconstructed. 

Key words: Methods: polarization - methods: data analysis - galaxies: magnetic fields - RM 
Synthesis - wavelets 



1 INTRODUCTION 

RM Synthesis is a new tool for the interpretation of polarized 
emission data in order to ge t information on the emitting media 
terentiens & de Bruvnll2005l) . Until now the method was mainly 
applied to di screte sources in Faraday space ("Faraday screens") 
(lHealdl l2009>). The aim of this paper is to discuss applications to 
sources which are extended in Faraday space. 

RM Synthesis is based on the fact that the complex polarized 
intensity P can be calculated as a Fourier transform 
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of the Faraday dispersion function F i n Faraday sp ace with the co- 
ordinate which is the Faraday depth (iBumll 19661) . Performing the 
inverse Fourier transform of P one obtains F which is the polarized 



intensity emerging from a region with Faraday depth 0. Faraday 
depth is defined as 

(p(x) = 0.81 I Bii(x')n^(x')dx\ (2) 
Jo 

where B\\ is the line-of- sight magnetic field component measured 
in jiG, Hq is the thermal electron density measured in cm~^ and the 
integral is taken from the observer to the current point x along the 
line of sight over the region which contains both magnetic fields 
and thermal electrons, where x' is measured in parsecs. A conven- 
tion is that is positive when B is pointing towards the Earth. 

In the context of RM Synthesis, one has to distinguish between 
Faraday depth and Faraday rotation measure RM which is defined 
as 

RMiA) = -g-, (3) 
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where ^ is the polarization angle (P = IPIe^^"^). 

The inversion of formula ([T) shows that the Faraday disper- 
sion function F is the Fourier transform of the complex polarized 
intensity: 



where k 



-p{k). 



(4) 
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20, and the Fourier transform is defined as 
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The practical limitation of the use of Eq. ^ comes from the fact 
that P can be measured only fo r /l^ > and in only in a finite 
spectral band (iBrown et al.ll2009h . 

The lower bound /Imin restricts the possibility to recognize ob- 
jects which are extended in Faraday space, the upper bound /Imax 
suppresses the visibility of small-scale structures of the object in 
Faraday space, and the lack of negative impedes the correct re- 
construction of the intrinsic polarization angle. Poor sampling in- 
side the observations window prevent the reconstruction of objects 
at large Faraday depths. All these problems and limitations can also 
be illustrated using the wavelet representation, which provides a 
kind of local Fourier transform, isolating a given structure on both 
sides of the Fourier transform. 

The topic of this paper is a further development of the wavelet 
approach as suggested by Frick et al. (2010) which presented at 
least a partial resolution of these difficulties. The paper is orga- 
nized as follows. In Sects. [2l3l we describe the wavelet-based RM 
Synthesis and the eff'ect of extending P to negative /l^, and in Sect|4] 
we compare RM-CLEAN (lHealdll2009h to a wavelet-based RM de- 
convolution. We subsequently describe a few simplistic, but rea- 
sonable models of magnetic fields and electron densities in galaxy 
discs. Large-scale features are described in Sect.|5] small-scale fea- 
tures in Sect.|7] The eff'ect of frequency coverage on reconstruction 
fidelity is discussed in Sect. [6l We subject these models to tradi- 
tional RM Synthesis and to our wavelet approach in Sect. [8] and 
summarize our results in Sect. |9] 



2 WAVELET-BASED RM SYNTHESIS 

The wavelet transform of the Faraday dispersion function F(0) can 
be written in the form 



(f}-b 



(6) 



where if/((p) is the analyzing wavelet, a defines the scale and b de- 
fines the position of the wavelet center. Then the coefficient Wf 
gives the contribution of the corresponding structure of scale b at 
position a to the function F. 

The contribution of the complex polarized intensity to the 
wavelet decomposition of the Faraday dispersion function can be 
divided in two parts wpia, b) - b) + w+{a, b), which are de- 
fined by P{X^) in negative and positive domain A^, correspondingly. 
From Eqs. ^ and ^ one gets 



oo 



(7) 



depth. Then, the coefficients b) can be reconstructed, follow- 
ing the symmetry arguments (see Eq. ([TSt below) 



w_(a, b) = w+ (a, 2(po(a, b) - b) , 



(8) 



where the function (po(a, b) takes on a value of b which corresponds 
to the maximum of the dominant structure near point (a, b). Distri- 
bution of \w+{a, b)\ is usually complex and applying of Eq. ^ can 
be nontrivial procedure. We suggest it in Secs.|6]and|8]how to pro- 
ceed in some particular cases. 

The function F can be restored using the inverse transform 



F(0): 



oo oo 



Wp (a, b) 
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(9) 



We use below the so-called Mexican hat 0(0) = (1-0^) exp(-0^/2) 
as the analyzing wavelet. The wavelet is real, however, the function 
P is complex, so the wavelet coefficients Wf are complex as well. 
For the chosen wavelet wpi-a, b) - wpia, b) and C^y - \. 



3 ROTATION MEASURE SPREAD FUNCTION 

The Faraday dispersion function F(0) obtained by RM Synthesis 
can be related to the true Faraday dispersion function F(0) by P(0), 
the "RM Spread Function" (RMSF) is 



F(0) = F(0)*P(0), 

where ^ denotes the convolution. The RMSF is defined as 



P(0) 



(10) 



(11) 



where W{A^) is the shape of the observable window in A^ space 
(the window function, which can include weights due to diff'erent 
sensitivities (gains) in the observation channels, see Heald (200i)) 
and ^ is a normalization constant. 

For a single spectral window (W = 1 for A^^^ < A < /Imax and 
W = elsewhere) the RMSF has the simple form 



P(0) = 



-2i(pAi 



sin (0A/12) 



(12) 



where /l^ 



('iii„+'^max)/2 and = (4,,-^^„)/2. R(4,) is shown 
in Fig. [it for A^ = 3.91 m^ and AA^ = 2.35 m^ (corresponding to 
the LOFAR wavelength range in the highband). Q The RMSF usu- 
ally is a complex- valued function which means that the obtained 
Faraday function diff'ers from the true one not only in amplitude, 
but also in phase. (Due the lack of data in the domain of negative 
A^ RM Synthesis rotate s the phase of the function under recon- 
struction.) To avoid this jBrentiens & de Bmvnl (l2005h proposed to 
multiply Eq. ^ by the factor e^'^^^o which means that 



/?55(0) = K J W(A^)Q-^''^^^'-^o^dA^ = K 



sin(0A/l2) 



and 



F(0) = e-2^^^5F(0)*P55(0). 



(13) 



(14) 



The coe fficients w+( a , b) ar e calculated from the known P(A^) for 
A^ > 0. iFrick et all bOlOh suggested to recognize the dominat- 
ing structures in the map \w+(a,b)\. The coordinate b of the corre- 
sponding maximum gives the position 0o of the structure in Faraday 



^ For the sake of definiteness, we consider here and below as a typical 
example the observational range in the LOFAR high band (1.25 < A < 2.5 
m, see Rottgering 2003) and use another observational range if required to 
illustrate particular properties of RM Synthesis. 
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Figure 1. RM spread functions for a single observation window with 1 .25 < 
A < 2.5 m. Panel (a): standard R; panel (b) Rbb', panel (c): Rw- Real part 
(thin solid red), imaginary part (dashed blue), modulus (thick black). 



Thus, the modified RM spread function Rbb((/>) has a vanishing 
imaginary part and a real part, which reproduces the shape of the 
envelope of the function R((/)) (Fig.lTJ)). 

Another way to improve the RM spread function follows from 
iFrick et aD fcOlOh . Their algorithm uses the symmetry argument 
(the results of Sect. 5 give strong support for this argument - real- 
istic objects mainly look like even objects in Faraday depth space). 
It means that if the object is centered at Faraday depth then 
F(20o - 0) = F((p). In space this gives 



(15) 



which means in terms of the spread function for a single window 

:)2^ 

(16) 



J? ^^ sin ((pAA^) . 
Rw((p) = cos (2(/)Aq) 



and 



F{cl>) = t-''^^'lF{cl>)^Rw((l>). 



(17) 



The RM spread function RwW is also shown in Fig.[T]:. Similar to 
Rbb this function has no imaginary part, but the real part remains 
exactly the same as in Eq. ([T2l) , compare thin solid lines in Fig[T^ 
and Fig[T]:. 

Note that the symmetry assumption used here works straight- 
forward only if there is only one source along the light of sight. If 
there are several sources the assumption should be applied using 
the wavelet technique. 



4 CLEANING AND WAVELETS 

We discussed above the reconstruction of objects symmetric in 
Faraday space around = 0o- The simplest objects for recognition 
are point-like sources, which can be described in Faraday space by 
delta functions 



F(0) = Foe2i^o^(0 - 0o). 

The corresponding polarized intensity is given by 



(18) 



(19) 



Note that the relation $T9\ determines P(A^) for A^ < and 
it follows the symmetry property ([T5] which is the starting point 
for our approach. It means that the traditional RM Synthesis for 
a point-like source indirectly implies this symmetry argument and 
can be considered as a particular case of the method presented in 
this paper. 

To recognize several point-like sources located on the same 
line of sight one can apply a deconvolution proced ure using th e 
Faraday dispersion function, called "RM-CLEAN" (lHealdll2009h . 
This procedure iteratively finds the maximum in the reconstructed 
Faraday dispersion function and then iteratively subtracts the scaled 
versions of the RMSF until the noise level is reached, after which 
a smoothed representation of the "CLEAN model" is used as the 
approximate true Faraday dispersion function. 

After removing the brightest structure (including the corre- 
sponding sidelobes) the procedure fin ds the maximum i n the re- 
maining data and restarts the iteration. iHeald et aP (l2009h showed 
how well RM-CLEAN deconvolved RM cubes from the WSRT 
SINGS survey. 

Let us show how this cleaning of objects, compact in Faraday 
depth space, can be realized in terms of wavelet-b ased RM Syn - 
thesis. We use the example, presented in Fig. 3from lH^(l2009h . 
This figure is partially reproduced in our Fig. |2] where the observed 
U(A^) and Q(A^), the restored F(0) before and after cleaning and 
the polarization angles observed and reconstructed after cleaning 
are shown. The observations are done in two windows. In Fig. [3^ 
we show the modulus of the wavelet coefficients calculated 
from same observed U (A^) and Q(A^). Then using the symmetry ar- 
gument for the main maximum we find and reconstruct F(0). 
Taking the dominating maximum in F(0) at = 0o = -185 rad 
m~^ we replace it by an equivalent point source and calculate for 
this point source the wavelet transform wi, based on the same sam- 
ple of observational channels as in original data. Then we subtract 
wi from the wavelet transform w of original data and we obtain the 
wavelet image, shown in the second panel of Fig.[3]3, which demon- 
strates a dominating maximum at = 165 rad m~^. Using again 
the symmetry argument we get w_ and restore the second point- 
like source in F ((/)). The corresponding polarization angles, calcu- 
lated for the given observational windows are shown in the panel 
(c), which is practically equivalent to the result of RM-CLEAN 
(Fig. [2]). Summarizing this section, we conclude that, in contrast 
to the extended sources, for a source that is point-like in Faraday 
space, wavelet-based RM Synthesis gives just the same result as 
conventional RM Synthesis. 



5 GALACTIC MAGNETIC FIELD IN FARADAY DEPTH 
SPACE 

The complex-valued intensity of polarized radio emission for a 
given wavelength is 
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Figure 2. Example of the RM-CLEAN process from lHealdl (l2009h for a 
bright point source in the field of the galaxy NGC 7331. Panel (a): ob- 
served Q (black) and U (green); panel (b): the restored F(0) before cleaning 
(green) and after cleaning (black); panel (c): observed polarization angles 
(green) and reconstructed after cleaning (black). 
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(20) 



defined by the emissivity s and the intrinsic polarization angle x 
along the line of sight. The integral is taken over the whole emit- 
ting region. The emissivity depends on the relativistic (cosmic-ray) 
electron density ric, the magnetic field component perpendicular to 
the line of sight and t he synchrotron spectral i ndex a (assuming 
a = 0.9 below) as fe.g. lGinzburg & Svrovatskiil(ll965h ) 



(21) 



Our analysis starts from an investigation of the typical shapes 
of the Faraday dispersion functions expected for lines of sight 
crossing galactic discs with several simple distributions of the par- 
allel magnetic field, which provides the Faraday rotation. These 
distributions are shown in Fig.l?^. We consider two symmetric dis- 
tributions with respect to center at x = xq. The first one has a sharp 
boundary at x = xq ± h. 



\x- 



xo\ ^ h 
xo\ > h 



(22) 



Here xq is position of galactic equatorial plane and h is the half 
thickness of the galactic disc. The second distribution has a smooth 
shape approximated by a Gaussian function with the scale height h. 
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Figure 3. Wavelet cleaning for the same example as in Fig. |2] Panel (a): 
wavelet plane after the first step (the symmetry of the main maximum is 
used); panel (b): wavelet plane after cleaning of the first point source and 
using symmetry argument for the rest (note that the amplitudes in the two 
figures are diff'erent); panel (c): observed polarization angles (green) and re- 
constructed after two-step wavelet cleaning (black). Here and below colors 
are given in a relative scale specific for each panel and black color corre- 
sponds to maximum. 



For comparison we take into consideration a distribution with 
a field reversal along the line of sight 



53||(x) = -as 



{X - Xq) 



-(x-xoY/h- 



(24) 



The amplitudes a2 and as are adjusted to obtain the same maximal 
Faraday depth for all three distributions. We choose an amplitude 
of the magnetic field of ai =2 yt/G and h = 0.5 kpc. The densities of 
relativistic electrons ric (responsible for synchrotron emission) and 
thermal electrons Hq (responsible for Faraday rotation) measured in 
cm~^ are assumed to have Gaussian profiles 



n^(x) = 0.03 e 



-(x-xoY/K 



(25) 



(23) 



Here C (measured in cm~^) is the relativistic electron density at the 
galactic midplane. We assume that the scale height of the relativis- 
tic electrons he is twice that of the thermal electrons and that 
the scale height of the magnetic field h is the same as that of the 
thermal electrons h^. We use C as a normalization constant for F 
and P, so that they are numerically evaluated in arbitrary but mu- 
tually consistent units. Fig. ^ shows that the symmetric profiles 
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Figure 4. Typical shapes of magnetic field distribution. Panel (a): 5i (solid 
line, red), B2 (dashed, blue), (dot-dashed, black); panel (b): correspond- 
ing Faraday depths 0. 
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Figure 5. Real part of the Faraday dispersion function F(0) (panel a), 
\mF = 0, and the modulus of corresponding PiA^) (panel b) for differ- 
ent magnetic field profiles: Bi (solid line, red), B2 (dashed, blue), B3 (dot- 
dashed, black). 



of the parallel magnetic field (first two cases) lead to a monotonic 
increase (decrease) of Faraday depth along the coordinate x inside 
the galaxy, while the antisymmetric profile (last case) leads to an 
inversion - the central galactic plane becomes the most distant (or 
most nearby) point in Faraday depth space. 

For the sake of normalization we choose the direction of the 
perpendicular field to obtain a purely real P (i.e. zero polarization 
angle). The perpendicular component is described by the same dis- 
tributions as for the parallel component in Eqs. (I22t-(l24l. Namely, 
we take Bi^ = ^ly, B2± - B2\\ and = B2\\. 
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Figure 6. Faraday dispersion function F(0) for different separations Ax be- 
tween the regions of synchrotron emission and Faraday rotation with Gaus- 
sian profiles: Ax = 0.1 kpc (solid line, red). Ax = 0.2 kpc (dashed, blue). 
Ax = 0.5 kpc (dot-dashed, black). 



Fig. |5^ shows that both even profiles {Bi and B2) produce in 
Faraday space basically box-like distributions with sharp bound- 
aries. In contrast, Faraday depth in the case of the magnetic field 
with a reversal varies slowly near the point where Faraday depth 
becomes minimal and has a sharp point-like structure in F(0) near 
the point where is maximal while F is still nonvanishing, in the 
given examples F is even maximal at this point. 

We conclude from Fig. that box-like and point-like shapes 
of F(0) are common in simple models of disc components of spiral 
galaxies. Of course, one can obtain a smoother distribution of F(0) 
by taking a wider distribution for thermal electrons than for rela- 
tivistic ones. According to Fig.[5]3 the level of polarized intensity at 
large is higher for pointed F profiles (hence easier to recognize 
by low-frequency observations) than for box-like ones. 

In general, the synchrotron emission is not necessarily located 
in the same region as the Faraday rotation. This leads to a so-called 
"Faraday screen" which can be due to a specific magnetic field con- 
figuration (strong weak B^) or to a difference between distribu- 
tions of nc{x) and nQ{x). For simplicity, we shift s{x) by some dis- 
tance Ax. The resulting Faraday dispersion functions for Gaussian 
shape (the case B2 in Fig. |4]) and different Faraday screen sepa- 
rations are shown in Fig. |6] For larger values of Ax the Faraday 
dispersion function F becomes more pointed and asymmetric. 



6 THE ROLE OF THE SPECTRAL RANGE 

The efficiency of RM Synthesis crucia lly depends on the obser - 
vational range /Imin < X < /Imax 

(iBrentiens & de Bruvnl l2005h . 
Wavelets give a helpful illustration for the role of these parameters. 

Fig. [7] shows the wavelet planes for the first artificial example 
(solid lines in Fig. [5]), calculated for different observational win- 
dows. Each panel shows the absolute value of the wavelet coeffi- 
cient W{a, b) in the (a, b) plane. The horizontal axis {b) presents 
the position in Faraday depth space, the vertical axis gives the 
scale a in logarithmic presentation. The upper bound is first fixed 
at /Imax - 2.5 m. The panel (a) shows the wavelet plane simulated 
for the lower bound /Lmin = 6 cm (thus 0.0036 < < 6.25 m^). 
The large horizontal dark feature at a ^ 40 rad m~^ corresponds 
to the box-like structure, and the two tapering structures are gen- 
erated by the sharp borders of this box. The next panel shows 
what remains in the wavelet plane if the lower bound is moved 
to /Imin = 21 cm (the window 0.044 < < 6.25 m^). The cen- 
tral spot (responsible for the box) is almost completely lost even 
at these relatively short wavelengths. Applying the actual LOFAR 
highband window (1.56 < < 6.25 m^), only weak traces of side 
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Figure 7. Wavelet plane \wF{a,b)\ of the example of Figs. IH and |5] (solid 
lines) for various spectral windows. Panel (a): 0.06 < /I < 2.5 m; panel (b): 
0.21 < i < 2.5 m; panel (c): \25 < A < 2.5 m; panel (d): 0.06 < A < 
0.21m. 



horns remain in the wavelet plane (panel c). The lower panel (d) 
in Fig. [7] shows the wavelet plane for a relatively high-frequency 
window (0.06 < A < 0.21 m), which perfectly keeps all informa- 
tion concerning the large scales, but completely loses any informa- 
tion concerning the small scales, responsible in this example for the 
abrupt boundaries. 

For applying the symmetry argument expressed by Eq. (|8) the 
position of each local maximum 0o and the domain of its influence 
should be defined. In Fig. [7^ we recognize three maxima and we 
divide the wavelet plane in three parts: one domain corresponds 




(a) 



(b) 



Figure 8. Test Faraday dispersion function (black dots in both panels), the 
same as shown in Fig. [5^ (solid red line), and results of RM Synthesis per- 
formed from the wavelet planes shown in Fig. |7] for various spectral win- 
dows. Panel (a): 0.06 < A < 2.5 m (solid, red line), 0.21 < A < 2.5 m 
(dashed, blue), 125 < A < 2.5 m (dash-dot, black); panel (b): 0.06 < A< 
1.0 m (solid, red line), 0.06 <A< 0.42 m (dashed, blue), 0.06 <A< 0.21 
m (dash-dot, black). 



a > 10 rad m~^ and two domains for a ^ 10 rad m~^ are separated 
at Z? ^ 18 rad m~^. Separatrices are shown in Fig. |7t by dashed 
lines. In the case shown in Fig. [7]: only two lower domains can 
be recognized. Note, that two horns-like structures are associated 
with the sharp borders, which are not symmetric objects, and the 
application of the symmetry assumption to them does not improve 
the reconstruction. 

Fig. [8] shows how the shrinking of the wavelength range of ob- 
servations reduces the quality of the Faraday function reconstructed 
from the corresponding sets of Wf(a,b). We see that increasing of 
/Imin leads to a gradual decay of the reconstructed F in the main 
bulk of the structure (Fig.[8t). Starting from /Imin = 21 cm we re- 
construct details associated with the boundaries only. Decreasing 
the upper bound /Imax we lose the small-scale details, restoring a 
smooth structure only (Fig.[8]D). 

The possibility to reconstruct a structure of scale A0 in Fara- 
day space is determined by the quantity X = A(pA^ where A is the 
wavelength at which the observations are performed. We illustrate 
the role of this quantity for the large-scale structure from Fig. [8^. 
The relevant parameters here are A0 = 38 rad m~^ and A - /Imin- 
X - 1.68 is for the dashed (blue) line and X = 60 for the dash- 
dot (black) lines. Because the latter curve keeps the information of 
sharp ends of the distribution F(0) and even the first curve looses a 
substantial part of the signal. 



X 



(26) 



is the condition to observe a galaxy as an extended source 
in Faraday space rather than as a couple of point-like sources 
terentiens & de Bruvnl l2005h . For observations in the LOFAR 
highband (/Imin = 1-25 m) this gives A0 = 0.64 rad m~^, which is at 
least one order of magnitude lower than a typical value for contem- 
porary galaxies. If /Imin = 0.167 m (by including observations e.g. 
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Figure 9. Modulus of the Faraday dispersion function F obtained from 
wavelet-based RM Synthesis performed for the observational window 
0.06 </l< 2.5 m with 1024 (panel a) and 128 (panel b) channels providing 
uniform sampling in space (dashed black line) or uniform sampling in 
frequency space (solid red line). 



with the ASKAP telescope) this yields a much more comfortable 
A0 - 35 rad m~^. 

For distant galaxies the wavelength in the framework of the 
source is /l/(l + z) (z is the redshift) and Eq. (|26] ) reads 



^ = A04J(1+z)2< 1. 



(27) 



Taking Z = 3 one obtains that A0 ^ 10 rad m~^ for LOFAR-type 
observations. LOFAR observations can also be used to isolate a 
contribution of weak intergalactic magnetic fields, which can give 
typical A0 much lower than that for galactic magnetic fields. 

As we conclude above, RM Synthesis of extended magnetic 
configurations requires a proper wavelength range of observations. 
/Imin has to be small enough to reproduce the grand design of F 
and a large /Imax is required to reproduce sharp details of F. Obvi- 
ously, RM Synthesis demands a substantial number of channels in 
the observational range to perform the Fourier transform included 
in the method. A dispersion function F which is extended over a 
large range of Faraday depths requires a large number of channels. 
If the number of channels is limited, the main question is how to 
distribute the channels between the given values of /Imin and /Imax- 
We clarify this using the example considered for Fig. |7] For the 
sake of the definiteness we compare two samplings, i.e. a uniform 
spacing in frequency space and that one in the space. 

We presume that the spectral range of observations is suffi- 
ciently large and use /Imin = 6 cm and /Imax = 2.5 m and apply 
the wavelet-based RM Synthesis. Fig. [9^ presents the results for 
a large (1024) number of channels. We see that even here a uni- 
form spacing in looks better than one in frequency space, which 
gives a more noisy result. The difference between the two kinds 
of spacings becomes much more pronounced for a moderate (128) 
number of channels. The results for the uniform spacing in is 
more or less similar to that one in the upper panel, while for a uni- 
form spacing in frequency the noise in F becomes comparable to 
the signal. It is clear that the eff'ect is essential for large spectral 
windows (namely, for a large ratio /Imax/^min) because a uniform 
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Figure 10. Magnetic field distributions over a galactic disc. The ratio of 
small-scale (b) and large-scale (B) components is b/B = 2, the ratio of 
spatial scales of the components is l/2h = 0.1. Panel (a): light-of-sight 
magnetic field; panel (b): polarized intensity. The long and short dashed 
lines give P for two different realizations of the turbulent magnetic fields 
calculated for a very narrow beam, while the solid line shows P calculated 
for a beam which contains 50 independent turbulent cells. 



spacing in frequency leads to points crowding at small A^. A poor 
sampling at large (essential to detect small scales in Faraday 
depth space) results in small-scale noise in the reconstructed sig- 
nal. Note th at the advantages of uniform spacing in A^ was first 
discussed bv lRuzmaikin & SokoloffI (Il979h . 



7 SMALL-SCALE MAGNETIC FIELDS 

In contrast to the traditional methods, wavelet-based RM Synthesis 
opens a new option to quantify the small-scale (turbulent) compo- 
nent of magnetic field. Remember that it is widely believed that 
one can isolate two contributions in total magnetic field, i.e. the 
large-scale (regular) magnetic field and the small-scale one. Cor- 
respondingly, the Faraday dispersion function F and the Faraday 
depth have two corresponding components, coming from large- 
scale and turbulent fields. The traditional Bum (1966) theory sug- 
gests that averaging P over turbulent variations yields an additional 
depolarization term like internal Faraday dispersion and beam de- 
polarization. In principle, one could try to modify Eq. ([T) to isolate 
the large-scale contribution and parameterize the small-scale one 
by intrinsic Faraday dispersion and beam depolarization. Wavelets 
allow us to perform this procedure in a more straightforward way 
and to obtain more detailed information concerning the small-scale 
component. 

We start here with an example of a magnetic field which con- 
tains large-scale and turbulent components (Fig. [TOh) and consider 
a slab with thickness 2h-\ kpc with a large-scale magnetic field B 
and a turbulent magnetic field with a Kolmogorov spectrum and an 
rms value b with the scale / = 0.1 kpc {ll2h - 1/20) and b/B - 2. 
The polarized intensity obtained from such a slab is aff'ected by 
two depolarization eff'ects. One eff'ect known as internal Faraday 
dispersion occurs because several (about 20 in our example) inde- 
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and various small-scale details at the bottom which represent the 
contributions of individual turbulent cells. Having a limited spec- 
tral range covered by the observations, we can obtain wavelet co- 
efficients in a limited part of the wavelet plane only. The border of 
this part is given by Eq. (|26] | which reads as 

Here we take into account that for structures of size a the range 
of Faraday depth is of order A0 = oT^ . The upper limit of the do- 
main admissible for LOFAR observations is shown in Fig. [TTIbv 
a horizontal dashed line. We conclude that the contribution of the 
large-scale magnetic field remains inaccessible for this type of ob- 
servations while the small-scale details remain visible. Note that 
the wavelet transform presented in the figure allows us to recog- 
nize the range of Faraday depths in which small-scale structures 
occur (from Orad m~^ to 38 rad m~^ for the given example). Obvi- 
ously, using data at large only, we cannot confirm (or reject) the 
existence of a large-scale magnetic field. 

The method under discussion presumes that the signal-to- 
noise ratio is sufficiently large. The lower panels of Fig.fTTb present 
the result of the wavelet RM Synthesis for P for which random in- 
strumental noise (independent in each spectral channel from the 
noise in the other channels) is added. The signal-to-noise ratio (cal- 
culated from the averaged signal and the noise per frequency chan- 
nel and per Q and U point for the observational range 1.25 < A < 
2.5 m) is about 2 (panel b) and the contribution of the small-scale 
fields is easily distinguishable from instrumental noise. These plots 
shows that a signal-to-noise ratio about 2 (per channel and per Q 
and U point) or higher is sufficient to isolate the range in Faraday 
space responsible for the small-scale magnetic fields. For illustra- 
tive purposes we give an example with the signal-to-noise ratio of 
about 0.5 and demonstrate that the contribution of small-scale field 
cannot be distinguished from instrumental noise (Fig.fTTt). 



Figure 11. Wavelet plane wpia, b) calculated for different ratios of signal- 
to-noise from P{A^), obtained for the wide beam (shown in Fig. 1 10b by 
the black line). From top to bottom: signal without noise, signal-to-noise 
ratios of 2 and 0.5. The horizontal dashed line shows the upper bound of the 
domain admissible to LOFAR-type observations. 

pendent turbulent cells are located at a given line of sight passing 
through the slab. This eff'ect occurs even for a very narrow beam, 
with a cross-section with the slab surface containing just one tur- 
bulent cell. The corresponding P for two independent realizations 
of the turbulent field is shown by dashed lines in Fig.fTOb. 

The other depolarization effect arises because the beam is usu- 
ally wide enough to contain many turbulent cells. If, say, the beam 
diameter on the slab is 1 kpc it contains about 100 independent 
cells. The polarized emission passing through independent cells 
contributes to P, leading to additional depolarization known as 
beam depolarization. The values of P calculated for a beam with 
a cross-section containing 50 independent cells is shown by the 
solid line in Fig.fTOb. We see that typical P/Pmax values are only a 
few percent for large A^. One could conclude that the emission is 
almost completely depolarized and useless for RM Synthesis. How- 
ever, the wavelet method being applied to P shows a well-ordered 
structure (Fig. [TT)- The point is that wavelets separate contribu- 
tions from different spatial scales and associate them with specific 
domains in the plot. 

In Fig. [TTt , we see a contribution from the large-scale mag- 
netic field (the long horizontal feature at the upper part of the panel) 



8 A SYNTHETIC EXAMPLE 

To summarize our findings and to conclude what can be learned 
from polarized radio observations of galactic magnetic fields using 
different RM Synthesis techniques, we now consider a complex test 
example. Suppose that we observe a galaxy which hosts both reg- 
ular and turbulent magnetic fields (like the example considered in 
Sect. 13 and a strong point-like source of polarized radio emission 
behind the galaxy on the same line of sight. Independent of its posi- 
tion in physical space (coordinate x) behind the galaxy, this source 
appears in Faraday space like attached to the galactic backside. For 
sake of definiteness we suppose that the contribution to F which 
comes from the galaxy is purely real and that one from the point 
source is purely imaginary. The task for the analysis is to recon- 
struct the Faraday dispersion model and then to restore the inten- 
sity and the polarization angle of the point source, to separate the 
contribution of large-scale and small-scale galactic magnetic field 
and recognize the shape of the large-scale structure. 

Fig.[T2k shows the wavelet plane calculated for this signal. We 
recognize here an extended structure in the region of large scales, 
10 ^ « ^ 50 rad m-^, which we identify with the contribution of a 
galactic disc as a whole, small-scale details at a < 1 rad m~^, which 
we identify as the contributions of the small-scale magnetic fields, 
and a bright horn located at ^ 37 rad m~^. Then we perform the 
wavelet-based RM Synthesis using the symmetry arguments to the 
galactic contribution with the reflection point at 0q - 18 rad m~^ 
and considering the domain of the horn as a point-like source, i.e. 
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Figure 12. Panel (a): Wavelet plane wpia, b) (in colors online) calculated 
from the polarized radio emission P{A^) generated by turbulent galactic 
magnetic fields superimposed on a point source behind the galaxy. The re- 
sults of RM Synthesis are shown in panels (b)-(d). Panel (b): wavelet recon- 
struction using the whole wavelet plane; panel (c): wavelet reconstruction 
using only the domains marked in the wavelet plane by dashed lines in panel 
(a); panel (d): reconstruction by standard RM Synthesis. The red thin line 
is for the real part, the blue thick line for the imaginary part, and the black 
dashed line for the modulus. 



apply the symmetry condition at the point 
means 



'po{a,b) = 



a> I 
a ^ I 



37 rad m ^. This 



(29) 



Applying the symmetry arguments for the whole domain would 
amplify the small-scale noise and reflect it with respect to the po- 
sition of the maximum. Therefore we perform a reflection locally 



(only in vicinity of the maxima) which leads to the definition 



W-(a, b) 



w+ (a, 2(po(a, b) ■ 
0, 
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(l)o(a, b)\ ^ a 
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(30) 



instead of Eq. The dashed lines in Fig.[T2b shows two domains 
around each maxima in which the symmetry arguments are used 
for calculation of w_ following Eqs. ( |29H30l ). 

The reconstruction in Fig.[T2t) shows a point-like source and a 
smooth contribution from the large-scale magnetic field, superim- 
posed onto fluctuations associated with the small-scale fields. If we 
are not interested in the small-scale details we ignore the wavelet 
coefficients outside isolated domains and obtain a smoothed F 
which represents the contribution of the mean field in the disc and 
the point- like source only (Fig.[T2b). We see a smoothed real contri- 
bution from the mean field and an imaginary point-like source, i.e. 
the method correctly reproduces the polarization angles. For com- 
parison with Fig. [T2b , we show in Fig. [T2h the result obtained by 
the standard RM Synthesis technique where the phase information 
is completely lost. More specifically, the blue thick line at Fig.[T2b 
follows the horizontal axis as required by the model while a nega- 
tive trend of the blue thick line is visible in Fig.[T2h. 



9 DISCUSSION AND CONCLUSIONS 

We considered how to apply RM Synthesis to extract information 
concerning magnetic fields of spiral galaxies from polarized emis- 
sion data. The main attention was given to two related novelties in 
the method, i.e. the symmetry argument and the wavelet technique. 
Symmetry arguments open the possibility to reconstruct P(A^) for 
negative from P(A^) for positive A^.lf F contains several contri- 
butions (say, the disc and an unresolved source) the reconstruction 
can be performed locally in Faraday depth space using wavelets. 
The wavelet technique appears to be useful in RM Synthesis in sev- 
eral areas which are not directly related to the problem of negative 
A^. In particular, they allow us to perform RM Synthesis locally in 
Faraday depth space. 

The idea of RM Synthesis was applied first for sources 
which are presumed to be point-like in Farad ay space or a 
combination of severa l point - like sources (iHaverkom et al .1120031: 
de Bmvn & BrentiensI l2005l : iHeald et all l2009l : I Wolleben et al. 



20101 : iLandecker et aLllioiob . We demonstrated that the traditional 
RM Synthesis for a point-like source indirectly implies a symme- 
try argument (see Eq. 22) and, in this sense, can be considered as a 
particular case of the method under discussion here. 

Investigating the applications of RM Synthesis to the polariza- 
tion details associated with small-sc ale magnetic fields, we isolated 
an o ption which was not c overed bv lBumI (ll966h . lBun] ([l966) (see 
also lSokoMetaP (Il998h ) described a contribution of small-scale 
fields in terms of Faraday dispersion and beam depolarization, i.e. 
using averages over a corresponding statistical ensemble. Here we 
use the complex polarization P for RM Synthesis without any aver- 
aging and demonstrate that it allows us to obtain ( if the A^ coverage 
is adequate) a range in Faraday space where the contribution from 
small-scale field is located. 

A general conclusion concerning the applicability of RM Syn- 
thesis to the interpretation of the radio polarization data for ex- 
tended sources like spiral galaxies is that quite severe requirements 
are needed to provide the full applicability of the method. The most 
severe one is the requirement fo r the wavelength ran ge where short 
wavelengths are highly desired. Ifirown et aD (l2009h show that RM 
Synthesis applied at the border of its applicability (say, if the num- 
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ber of frequency channels is very small) can sometimes be mislead- 
ing. 

Traditional methods of the pattern recognitio n for mag - 



netic structure s in spiral g alaxies (e.g. iRuzm aiki n et al 
Sokoloff etal.1 (Il992): Fric k et all (j^OOOh 



Fletcher et all (120041) : 



1990); 



Patrikeev et^ (I2OO6) : iStepanov et al.l (l2008h ) were mainly based 



on the interpretation of polarization angles at a few wavelengths 
only. The efforts of the above papers concentrated mainly on the 
distribution of polarization angles in the radio image. Modern wide- 
band observations with many frequency channels and RM Synthe- 
sis open wide perspectives to extract much more information from 
the polarization signal in a given telescope beam. However, they 
cannot fully replace the traditional analysis of the distribution of 
the polarization patterns. Of course, the development of RM Syn- 
thesis in order to deal with radio images as a whole rather than with 
just one beam is a goal for further developments, and wavelets are 
likely to be helpful here. However, this is outside of the scope of 
this paper. 

Obviously, a quantified characteristic of reconstruction fidelity 
of the algorithm under discussion is highly desirable. However, a 
general norm is not informative in a spectral problems with varying 
windows - one can completely lose the large-scale structures and 
reconstruct precisely the small-scale details (or vise versa). Fig. [8] 
is a good illustration for that: the central part of the structure is 
completely lost (for 0.06 < /I < 0.42), but the borders are recon- 
structed quite well. Thus the "global fidelity" is very poor, while 
the "fidelity of border reconstruction" is fine. So, we avoided inten- 
tionally to quantify the quality of the reconstruction by a norm in 
some functional space. 

The RM Synthesis approach was suggested in the context of 
observations with a rather large value of /Imin while the traditional 
method to measure RMs between only two wide frequency bands 
intensively used the data at /I = 3 cm and X = 6cm and rarely used 
the wavelength range X > 22 cm. The attention of the traditional 
method was focused on the magnetic field configuration rather than 
on the Faraday dispersion function which appears in this context as 
an intermediate rather than a final result. On the other hand, ob- 
servations in the wavelength range 1.25 < A < 2.5 m (the LOFAR 
high band) cannot recognize a lot of structures which would be ac- 
cessible to observations in the wavelength range 0.03 < X < 2.5 m 
with an adequate coverage. However, this option will have to wait 
for future radio telescopes like the SKA. 

This paper, although based on model examples, presents an 
optimistic view that the combination the ideas of traditional RM 
studies of galactic magnetic fields at short wavelengths and those of 
RM Synthesis at large wavelengths is possible. The final conclusion 
can be made after application to real data. 
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